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ABSTRACT 

Using the nonperturbative Schwinger-Dyson equation, we show that chi- 
ral symmetry in weak-coupling massless QED is dynamically broken by a 
constant but arbitrarily strong external magnetic field. 
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Chiral symmetry plays an important role in elementary particle and nu- 
clear physics. In this Letter we examine its breaking in the theory of quan- 
tum electrodynamics. It has been known 0] for some time that QED may 
have a nonperturbative strong-coupling phase, characterized by spontaneous 
chiral symmetry breaking, in addition to the familiar weak-coupling phase. 
The existence of this new phase was exploited in a novel interpretation of 
the multiple correlated and narrow-peak structures in electron and positron 
spectra||| observed at GSI several years ago. According to this scenario, the 
e + e~ peaks are due to the decay of a bound e + e~ system formed in the new 
QED phase, which is induced by the strong and rapidly varying electromag- 
netic fields present in the neighborhood of the colliding heavy ions. While 
the experimental situation with regard to these anomalous e + e~ events is un- 
clear, especially after similar experiments at Argonne have yielded negative 
resultsfU, it is still of great interest to investigate whether background fields 
can be physically used to induce chiral symmetry breaking. Now the question 
is: what kind of background fields can potentially induce chiral symmetry 
breaking in gauge theories? We recall that in a magnetic monopole field a 
gauge field breaks chiral symmetry and that in the Nambu-Jona-Lasinio 
model a magnetic field drives the critical transition point towards weaker 
couplingf?]]. Thus, magnetic fields are obvious candidates. We will take ad- 
vantage of the strong-field techniques introduced by Schwinger and others 
to consider a constant magnetic field of arbitrary strength. In order to put 
our problem in as general a setting as possible, we will use the nonperturba- 
tive Schwinger- Dyson equation approach. We will then compare our result 
with that obtained recently by Gusynin, Miransky, and Shovkovy|8|] whose 
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method is very different from ours. A comparison of the two approaches 
will sharpen our understanding of the underlying physics and the kind of 
approximations involved. 

The motion of a massless fermion of charge e in an external electromag- 
netic field is described by the Green's function that satisfies the modified 
Dirac equation proposed by Schwinger: 

7 • Il{x)G A (x, y)+ J d 4 x'M{x, x')G A (x', y) = 6® (x-y), (1) 

where H^x) = —id^ — eA^x), and M(x,x') is the mass operator M in 
the coordinate representation. For a constant magnetic field of strength H, 
we may take A 2 = Hx\ to be the only nonzero component of A^. In the 
following we will use the method due to Ritus0, which is based on the 
use of the eigenfunctions of the mass operator and the diagonalization of 
the latter. As shown by Ritus, M is diagonal in the representation of the 
eigenfunctions E p (x) of the operator (7 ■ II) 2 : 

-( 1 -U) 2 E p (x)=p 2 E p (x). (2) 

The advantage of using this representation is obvious: M can now be put 
in terms of its eigenvalues, so the problems arising from its dependence on 
the operator II can be avoided. In the chiral representation in which 03 
and 75 are diagonal with eigenvalues a = ±1 and x = ±1, respectively, the 
eigenfunctions E pax (x) take the form 

E pax (x) = Ne^+^+^DMu^ = E pax u ax , (3) 

where D n (p) are the parabolic cylinder functions [0 with indices 

eHa 1 

n = n(k, a) = k + —— - -, k = 0, 1, 2, (4) 
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and argument p = \j2\eH\(x\ — ^). Note that n = 0, 1, 2, ... . The nor- 
malization factor is N = (ATt\eH\) l l A / \fn}.\ p stands for the set (po,P2,P3, k); 
and u) ax are the bispinors of 03 and 75. 

Following Ritus, we form the orthonormal and complete []TT|] eigenfunction- 
matrices E p = diag(-E p ii, E p _ u , E pl _ ly They satisfy 

7-n E p (x) = E p {x) 7 .p (5) 

and 

M(x, x')E p (x') = E p (x)5 {A) (x-x')t A (p), (6) 

where represents the eigenvalues of the mass operator, and po = 



Po, pi = 0, p2 = —sgn(eH)^j2\eH\k, p 3 = p%. These properties of the 
E p (x) allow us to express the Green's function and the mass operator in the 
.^-representation as (E p = 7°_E^7°) 

G A (x, y) = f -^-E p (x) _ X ~ E p (y), id A p = J2 I dp dp 2 dp 3 , 

(?) 

and 

M(p,p) = [ d A xd A x'E p (x)M(x,x')E p ,(x') = t A (p){2-n) A 5 (A \p - p), (8) 



respectively, where 5^ A) (p - p') = 5 kk >5(p - p' )S(p 2 - P 2 )S(p 3 - p' 3 )- 
We work in the ladder approximation in which 

M(x, x r ) = ie 2 YG A (x, x'^D^x - x'), (9) 

where D^ u (x — x') is the bare photon propagator, 



The Schwinger- Dyson (SD) equation then takes the form 

t A {p){2^5{p-p') = ie 2 j ' d A xd^x'f^-E p {x)^E^{x) 

' ^-E^ix'WE^D^x-x'). (11) 



7-p"+E A (p' : 

After integrations over x, x', p' ', p' 2 ', and p' 3 ', the SD equation is simplified 



to read (r = J (qf + q%)/(2\eH\), <p = tan^^/q^) 



d A q 1 



X A (p)5 w = ze 2 ^ / 7 ^^====e r e 

k" 



r sgn(e_ff)(n'— n+n"— n")i/3 



(2tt) 4 Vnln^n'Mn"! 



x1(^-(i-0^)t°AtVa" 

x l -~ 7° A" 7 7 i/ A'J nn " (r) J fi . B , (r) , (12) 

7-p + X A (p ) 

where summing over cr, a', a", and a" on the right hand side is understood, 
and 

min(n,n') | /| 

Jnn>(r)= — ^J- ^sgn(eH)rT +n '- 2m . (13) 

^ ml (n — my. (n' — m) ! 



We have also used the following notations|12| in Eq. (|I2D : p" = p — qo, p" = 0, 
p" = - sgn(eif) J2|eif|/c", p" = p 3 -<?3, A = A(cr) = diag(5 (T i, ^-i, 5^, 5 CT _i), 



A' = A(tr'), - etc., ri = n{k',a'), n" = ra(fc", a"), and n" = n(Jfe",cr"). 

Eq.(^) may be solved by following the standard procedure flBf of writing 



Ea(p) = /3j ■ p + Xa(p)> where S^(p) corresponds to the dynamically gen- 
erated fermion mass. We will assume that £a(p) is proportional to the unit 
matrix (it will be seen later from the solution that this is a self-consistent 
assumption). Eq. ([T2|) then leads to two coupled equations for (3 and S^: 

\ 8 = i£ 2y f A 1 p -ry sgn( e H)(n'-n+n»-n") y 

^ ^ 7 .p J (27r)Vn!n'!n"!n"! 




(1 + /3)V 2 + ^ 
f { -(1 + /?)(G 2 -^Q 2 ) ) 



(14) 



where G x = A7^A"A" 7At A' = -2(5 a * 1 5s»i+5 »-. 1 5 s »- 1 )difLg(5 al 5 a , 1 , 

<ki<W <J ff -i<f«*-i), Qi = A( T -g)A"A"( 7 -g)A / , G 2 = A 7 "A" (7 -p^A^A', 

and Q 2 = A( 7 • g) A" (7 • p 5 ' ) A" (7 • q) A'. 



We seek solutions with (3 — 0. We will show later that such a solution 
is consistent only with the Feynman gauge (£ = 1). In this case the two SD 
equations decouple and only G\ is relevant for determining the dynamical 
fermion mass. The spin structure of G\ implies that necessarily cr" = <t", 
which, in turn, implies that necessarily n" = n" . It is convenient to make a 
change of integration variables from (q\, q 2 ) to the "polar coordinates" (r, tp). 
The integration over (p yields 



We note that the spin structure of G\ also implies that a = a', which, 
together with the 5 nn i from Eq. (|15|) , matches the 5kw on the left hand side 



Due to the factor e _r in the integrand in Eq.([L4]), contributions from 
large values of r are suppressed. Let us therefore, as an approximation (we 
will find out later what physical condition validates this approximation), keep 
only the smallest power of r in J nn »(r), i.e., 




(15) 



of Eq.flTJ). 



J nn" \J) 



[max(n, n")]! 



(i sgn(e#)r)l n - n "L 



(16) 
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Since the leading contributions come from the term corresponding to n" = n, 
we need only keep the term given by k" = n+\ — ^-sgn(eif ) in the summation 
over k" . As a result, we can replace J nn » by n\. The SD equation (Eq.flTJ)), 
thereby vastly simplified, becomes 

*m - ^h^r^-'^^h (i7) 

where q 2 = —q 2 . + gf + 2\eH\r 2 and p" 2 = — (p — go) 2 + (P3 — qz) 2 + 2\eH\k" . 

Let us make a Wick rotation to Euclidean space: po — > ip^, go — > i^i- 
Consider the case with p = 0, i.e., po = p^ = k = 0. Notice that k = means 
that, for positive (negative) ei/, a = 1(— 1) and n = 0, the last of which 
implies that k" = and cr" = 1(— 1) for the respective sign of eH. We also 
note that for either sign of eH, the matrix G\ can be effectively replaced by 
—2 x 1. We will assume that the dominant contributions to the integral in 
Eq.flTT|) come from the infrared region of small q% and q^ (this assumption 
will be seen to be self-consistent). Thus, it is reasonable to replace XU(p") 
in the integrand by S^(0) = m x 1. Eq.([T7|) then becomes 

m ~ ~2 I d( l3dq4 I dr 2 e~ r \ 1 2 2 ™ 2 (18) 

7r 2 J Jo 2r 2 + L 2 (q$ + qf) m 2 + (gf + q%) 



where a = e 2 /47r is the fine structure constant and L = l/^J\eH\ is the 
magnetic length. The integrations over q^ and q^ give 

a r°° 2 e~ r2 \n(2r 2 / m 2 L 2 ^ 



!--/ dr O 2 2T2 ( 19 ) 

n Jo 2r z — m z L z 



which yields the nonzero dynamical mass as 



m ~ a 



eH\ e _ VlF^ (20) 



where a, b, and c are constants of order 1. 

Eq.(^) clearly demonstrates the nonperturbative nature of the result. 
It also shows that our approximations break down when a > 0(1). As a 
further check on the consistency of our assumptions, we note that, according 
to Eq . (|18|) , the dominant contributions to the integrals come from the region 
2r 2 ~ m 2 L 2 ~ L 2 (q 2 + Our earlier assumption that effectively r <C 1 
is now translated to the physical assumption that mL <C 1, which requires 
that a <C 0(1); in other words, the dynamical chiral symmetry breaking 
solution we have found applies to the weak-coupling regime of QED! Now it 
is also evident that indeed the infrared region of q^ and q^ gives the dominant 
contributions to the integrals. 



It remains for us to show that (3 = solves Eq. (|14]) only if £ = 1. The 
main point to note is that, consistent with the k" = approximation made 
above, we can approximate the j-p" in G2 by — (7°go + 1 3 <13) (recall also that 
we are considering the case of p = 0). But then the piece of the integrand 
involving G2 is odd in q as well as in g 3 , and hence vanishes upon integration. 
It follows that the solution (3 = requires £ = 1, i.e., the Feynman gauge. As 



a result, the Qi-piece on the right hand side of Eq. (14) does not contribute. 
(That is fortunate because Q\ is actually not proportional to the unit matrix; 
its presence in the SD equation would have spoiled the assumption that 
is the dynamical mass multiplied by the unit matrix.) 

In summary, we have found a solution to the Schwinger- Dyson equation in 
the presence of an arbitrarily strong constant magnetic field, which indicates 
that, even at weak gauge coupling, an external magnetic field can trigger 
the dynamical breaking of chiral symmetry in QED, with the dynamical 



mass of the fermion given by Eq. (p0|) . Our general conclusion agrees with 
a recent finding by Gusynin et a/.||, whose approach is very different from 
ours. It would be interesting to examine if there are additional solutions 
of chiral symmetry breaking due to an external magnetic field. A parallel 



calculation for the case of a constant background electric field [14j] or other 
background field configurations may also shed light on the dynamics of chiral 
symmetry breaking in gauge theories. The formalism proposed here will be 
most suitable for these studies. 
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